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We study and extend the semidefinite programming (SDP) hierarchies introduced in [Phys. Rev. 

Lett. 115, 020501] for the characterization of the statistical correlations arising from finite dimen¬ 
sional quantum systems. First, we introduce the dimension-constrained noncommutative polynomial 
optimization (NPO) paradigm, where a number of polynomial inequalities are defined and optimiza¬ 
tion is conducted over all feasible operator representations of bounded dimensionality. Important 
problems in device independent and semi-device independent quantum information science can be 
formulated (or almost formulated) in this framework. We present effective SDP hierarchies to at¬ 
tack the general dimension-constrained NPO problem (and related ones) and prove their asymptotic 
convergence. To illustrate the power of these relaxations, we use them to derive new dimension wit¬ 
nesses for temporal and Bell-type correlation scenarios, and also to bound the probability of success 
of quantum random access codes. 

PACS numbers: 


I. INTRODUCTION 

Many problems in quantum information theory can be 
formulated as optimizations over operator algebras of a 
given dimensionality. Let us quickly review some of them. 

In 1-way quantum communication complexity [1, 2], 
two separate parties, call them Alice and Bob, are re¬ 
spectively handed the bit strings x, y £ {0,1}". Bob’s 
task consists in guessing the value of the Boolean func¬ 
tion f{x,y) £ {0,1}, and, to this aim, we allow Alice 
to send him a D-dimensional quantum system. Under 
these conditions, computing the maximum probability 
that Bob’s guess is correct amounts to optimizing over 
all possible D-dimensional quantum states prepared by 
Alice and over all possible measurements conducted by 
Bob on such states. 

In Bell scenarios [3], two or more distant parties con¬ 
duct measurements over an unknown quantum state. It 
has been observed that, even if we do not assume any 
knowledge whatsoever about the mechanisms of the mea¬ 
surement devices, it is sometimes possible to lower bound 
the dimensionality D of the quantum systems accessible 
to each party by virtue of the correlations between the 
measurement results alone [4-7] . In this regard, deriving 
dimension witnesses , i.e., statistical inequalities satisfied 
by the correlations achievable through multipartite quan¬ 
tum systems of local dimension D, can be understood as 
an optimization over entangled states and measurement 
operators. 

Entanglement distillation [8], or the capacity to pre¬ 
pare states close to a pure singlet given a number of 
mixed states through Local Operations and Classical 
Communication (LOCC), is one of the most conventional 
problems in quantum information science. More gener¬ 
ally, determining whether the state transformation p —> a 
can be effected via LOCC can be interpreted as a feasibil¬ 
ity problem, where the free variable is the corresponding 


LOCC map. If we restrict to local protocols or 1-way 
LOCC, the set of relevant maps admits a simple charac¬ 
terization in terms of tensor products of Krauss operators 
satisfying certain quadratic constraints. 

The above problems involve optimizations over a tuple 
of noncommuting variables Xi,...,A ra satisfying a num¬ 
ber of polynomial constraints, such as Xf = X t = X- (for 
projectors), or XiXj = X\X t = I (for unitaries). The 
number of total constraints is typically so low that, even 
fixing the dimensionality D of the spaces where these 
operators act, we find a continuum of inequivalent repre¬ 
sentations. 

Analogous problems emerge in the black-box approach 
to quantum information theory [3, 9-11], where the only 
constraints considered are essentially commutation rela¬ 
tions between projection operators implemented by dis¬ 
tant parties. The characterization of quantum nonlocal¬ 
ity has boosted the field of noncommutative polynomial 
optimization (NPO) theory [12-14], where the goal is, 
precisely, to conduct optimizations over all tuples of op¬ 
erators satisfying a number of polynomial inequalities. 
NPO theory achieves this via hierarchies of semi-definite 
programming (SDP) relaxations whose first levels ap¬ 
proximate quite well the space of feasible solutions. 

Unfortunately, NPO theory does not offer any means 
to bound or fix the dimension of the Hilbert spaces where 
such operators act. Since the afore-mentioned problems 
in quantum information theory become senseless or triv¬ 
ial in the high dimensionality limit, one would not expect 
NPO to be of any use. 

This view changed with the publication of Ref. [15], 
where a systematic way to devise hierarchies of SDP re¬ 
laxations for a wide class of NPO problems under dimen¬ 
sion constraints was introduced. Such relaxations, which 
seem to work quite well in practice, were used to de¬ 
rive a number of new results in quantum nonlocality and 
quantum communication complexity. Important theoret- 
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ical aspects, such as the completeness of the hierarchies, 
or the explicit nature of the dimension constraints, were 
nonetheless left out. Actually, from a reading of Ref. [15], 
it is not even clear which problems can be attacked with 
the new tools. 

In this paper, we generalize the SDP schemes proposed 
in Ref. [15] to cover all NPO problems where the di¬ 
mensionality of the relevant Hilbert spaces is bounded. 
We will prove the convergence of the resulting SDP hi¬ 
erarchies and discuss their efficient implementation. Fi¬ 
nally, we will use them to derive a number of new re¬ 
sults in quantum information theory: from new bounds 
for quantum random access codes (QRAC) for both real 
and complex quantum systems to semi-device indepen¬ 
dent Positive Operator Valued Measure (POVM) detec¬ 
tion [16, 17]; from the characterization of temporal cor¬ 
relations [18] under dimension constraints to the explo¬ 
ration of tripartite Bell scenarios where the dimensional¬ 
ity of just one of the parties is limited. 

The structure of this paper is as follows: first, we will 
define the generic problem of NPO under polynomial con¬ 
straints. Then, in Section III, we will present a hierar¬ 
chy of SDP relaxations to tackle it. In Section IV we 
will prove the convergence of this hierarchy. As it turns 
out, a straightforward implementation of the hierarchy 
would converge too slowly to be of much use, given rea¬ 
sonable computational resources. Hence in Section V we 
will give some hints to boost the speed of convergence 
of the method—to make it practical--, and exemplify 
its application by solving a specific problem on temporal 
correlations [18]. In Section VI we will explore the perfor¬ 
mance of related SDP hierarchies to characterize quan¬ 
tum nonlocality under dimension constraints and quan¬ 
tum communication complexity. In Section VII we will 
offer some advice on how to code the corresponding pro¬ 
grams. Finally, we will summarize our conclusions. 


II. NONCOMMUTATIVE POLYNOMIAL 
OPTIMIZATION UNDER DIMENSION 
CONSTRAINTS 


where the maximization is supposed to take place over 
all Hilbert spaces hi with dim(H) < D , all tuples of op¬ 
erators (X\, ...,X n ) C B(H) and all normalized states 
l^) 6 hi. Note that, if it were not for the dimension 
restriction D, the above would be a regular noncommu- 
tative polynomial optimization problem [14]. 

We say that the relations 1Z satisfy the Archimedean 
condition if there exist polynomials fj(X), gij(X) such 
that 

* 3 

+ Y / 9ij(X)Ui(X)g ij (X). ( 2 ) 

In the following, we provide a hierarchy of semidef- 
inite programming relaxations for this problem. Such 
a hierarchy will provide a decreasing sequence of values 
Pi > P 2 > such that p k > p*,\/k. Moreover, if the 
Archimedean condition holds [68] , then the hierarchy can 
be shown complete, i.e., limfe_ >00 p fc =p *. 

III. THE METHOD 

Let y = (y w )\ w \< 2 k be a sequence of complex numbers 
labeled by monomials w of the variables Xi,...,X n of 
degree |ui| smaller than or equal to 2k. Such a sequence 
will be called a 2 k th -order moment vector. Given y , the 
k th - order moment matrix M k (y) is an array whose rows 
and columns are labeled by monomials of Xi,...,X n of 
degree at most fc, and such that 

Adk{y)u,v = Vul v ( 3 ) 

Given y = (y w )\w\< 2 k and a Hermitian polynomial 
q(X) = QwUi(X), where the w in the summation 
ranges over all monomials of Xi, ..., X n of degree at most 
deg(g), the corresponding fc -order localizing matrix is 
defined as 


Consider the set S of all n-tuples of self-adjoint oper¬ 
ators (Xi,... ,X n ) satisfying the relations 1Z = {qi(X) > 
0:i = l,..., to}. Here qi(X) denotes a Hermitian polyno¬ 
mial of the variables Xi,.... X n , while the notation A > 0 
signifies that operator A is positive semidefmite. We will 
call each feasible tuple (Xi,..., X n ) a representation of 
the polynomial relations 1Z. Given a Hermitian polyno¬ 
mial p(X) and a natural number D, the problem we want 
to address is how to maximize the maximum eigenvalue 
of p(X) over all representations of 1Z of dimension D or 
smaller. In other words, we want to solve the problem 


p* = max (ib\p(X)\xb) 

s.t. dirn(TI) < D,qi{X) > 0, for i = 1, ...,to, (1) 


M k (qy) u ,v — 'y \ qwPu\u 


( 4 ) 


with |u|, |u| < k — 

A sequence y = (y w )\ w \< 2 k admits a quantum repre¬ 
sentation if there exists a representation (X\, ...,X n ) C 
B(H) of relations { qt(X ) > 0}i, with dim(TI) < D , and a 
normalized vector \ijd) S hi such that y w = {'i/j\w(X)\'i/j). 
It is a standard result in non-commutative polynomial 
optimization theory that, if (yw) admits a moment repre¬ 
sentation (of whatever dimensionality), then M k (y) and 
M k (qiy) must be positive semidefmite matrices for all 
orders k [14], 

The above positive semidefmite constraints are not 
dimension-dependent, and are actually obeyed by mo¬ 
menta emerging from representations of {qi(X) > 0}^ 




3 


of arbitrary (even infinite) dimensionality. The key to 
introduce dimension constraints is to acknowledge that 
moment vectors (y w ) admitting a quantum representa¬ 
tion satisfy a number of extra linear restrictions depend¬ 
ing on the value of D. 

Some of such restrictions arise due to matrix poly¬ 
nomial identities (MPI) [19]: these are polynomials 
s(X) of the variables Xi,...,X n which are identically 
zero when evaluated on matrices of dimensionality D 
or smaller. For D = 1, all MPIs reduce to commuta¬ 
tors, i.e., [Ai,Aj\ = 0, if Ai,Aj G B( C). Identifying 
At = Xi, this implies that sequences y = (y w )\ w \< 2 k 
admitting a 1-dimensional moment representation must 
satisfy yx x x 2 ~ Vx 2 x x = 0. Actually, for any value of D 
there exist MPIs from which non-trivial linear constraints 
on y can be derived. For D = 2, all MPIs are gener¬ 
ated by composition of the identities [[Ai, A 2 ] 2 , A 3 ] = 0 
and E 7r e5 4 s g n ( 7r )^(i)^(2)^(3)^(4) = 0, where S 4 
denotes the set of all permutations of four elements. The 
latter identity is a particular case of the family of poly¬ 
nomial identities Id, with 

^2 sg n ( 7r )^(i)---A 7r ( d ) =0. (5) 

eSd 

It can be proven that all D x D matrices satisfy I 2 d [19], 
also called the standard identity. The problem of de¬ 
termining the generators of all MPIs for dimensions D 
greater than two is, however, open. 

A non-trivial relaxation of problem (1) is thus: 

p k =ma p w y w , 

y z ' 

w 

s.t. y G Sd,W = 1, M k (y) > 0, 

M{qiy)> 0, for * = 1, (6) 

where Sp denotes the span of the set of feasible sequences 
y = {yw)\w\< 2 k- This is a semidefinite program, and, 
as such, can be solved efficiently for moment matrices 
of moderate size (around 200 x 200) in a normal desk¬ 
top [20]. 

Equivalently, we can re-express the positivity condi¬ 
tions as M/c = M k {y) ® ®™ 1 Mfc(g i y) > 0 and rewrite 
the objective function as a linear combination of the en¬ 
tries of the first diagonal block of M. That way, we can 
regard the block-diagonal matrix M (and not y) as our 
free variable, hence arriving at the program 

p k = max y - ' p w M w i, 

s.t. Me M k D ,Mi,i= 1,M> 0, (7) 

where M ^ denotes the span of the set of feasible extended 
moment matrices. This reformulation of problem (6), al¬ 
though conceptually more cumbersome, leads to simpler 
computer codes. 


The key to implement either program is, of course, to 
identify the subspaces 5^,, M k D . We will now provide 
two methods to do so. Both have advantages and disad¬ 
vantages. In [21] we provide yet a third method, which, 
although more complicated than the other two, requires 
considerably less memory and time resources, making it 
suitable for high order relaxations. 

The randomized method 

We sequentially generate n-tuples of random Her- 
mitian D x D complex matrices X 5 = (X(, ...,X^) 
and normalized random vectors \ip^) G C D , which 
we use to build moment and localizing matri¬ 
ces M{ v = {^\u{X^v{Xi) |V®, M( qi )i tV = 
(ipi \u{X^yq i {X)v(X : >)\ip^), respectively. Their direct 
sum will constitute an extended moment matrix M 3 . 
Adopting the Hilbert-Schmidt scalar product ( A, B) = 
tr(Ad.B), one can apply the Gram-Schmidt process [69] to 
the resulting sequence of feasible extended moment ma¬ 
trices in order to obtain an orthogonal basis M 1 , M 2 ,... 
for the space spanned by such matrices. We will no¬ 
tice that, for some number N , M N+1 = 0, up to nu¬ 
merical precision. This is the point to terminate the 
Gram-Schmidt process and define the normalized matri¬ 
ces {r* = ' 3 = It ea sy to see that, 

even though the matrix basis } was obtained ran¬ 
domly, the space it represents is always the same, namely, 
M k D . 

Indeed, let N = dim(Af)j), and suppose that 
M 1 ,..., M- 7-1 are non-zero, with j < N. Then the en¬ 
tries of the matrix M J will be polynomials of the com¬ 
ponents of X^i/ji. Since j < N, there exists a choice 
of 2 s such that M J (X J , -i/d) ^ 0. The probability that a 
non-zero polynomial vanishes when evaluated randomly 
is zero, and so we conclude that M° will be non-zero 
with probability 1. On the other hand, M 1 = M ^ 0, so 
by induction we have that N randomly chosen moment 
matrices will span Ad k D with certainty. Consequently, 
M n + i = 0 indicates when to stop the procedure. 

Remark 1. A cautionary note is in order: for high 
orders k, it is expected that program (7), as it is writ¬ 
ten, will not admit strongly feasible points. That is, the 
subspace M k D will not contain any positive definite ma¬ 
trix. This can be problematic, as many SDP solvers need 
strong feasibility to operate. The solution is to add up 
the random extended moment matrices M J as we produce 
them, i.e., to compute the operator T = ■ 

Since { M J } were randomly generated, it can be argued 
that the support of any matrix M G is contained in 
the support ofT. Let V be any matrix mapping supp(T) 
to C dlm ( supp ( Td isometrically. We just need to replace 
the positivity condition M > 0 in (7) by VMV t > 0, 
which, by definition, admits a stricly feasible point. 

This method has the advantage that it is extremely 
easy to program; more so when the constraints {qfiX) > 
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0} reduce to polynomial identities, as we will see. One 
disadvantage is that, in practice, the decision to stop 
the protocol amounts to verifying that the entries M N+1 
are zero up to e precision. Choosing the right value for 
the threshold e is a delicate matter: if too small, the 
algorithm will not halt; if too large, the algorithm will 
stop before it finds a complete basis for M k D . The second 
problem is that, due to rounding errors, it is possible that 
the algorithm will not identify the right subspace, but 
only an approximation to it. 

The deterministic method 

We choose a simple distribution f(X,ip)dXdip, say, 
a Gaussian, for the entries of each of the matrices 
X \,..., X n and the components of the unnormalized vec¬ 
tor ip. Then we define the components of the 2k moment 
vector y(X,ip) via the relation y(X,ip)w = {ip\w(X)\ip) ■ 
Then we compute analytically the matrix 

s = J f(x, <p)dXdMx, ip)y(x, ( 8 ) 

Clearly, the space S|> corresponds to the support of S. 
Diagonalizing S and keeping just the eigenvectors with 
non-zero eigenvalue we hence obtain an orthonormal ba¬ 
sis for Sp. 

The disadvantage of this method is that it involves 
symbolic computations, and hence, depending on the 
platform used, it is either more difficult to code or re¬ 
sults in slower programs. 

IV. CONVERGENCE OF THE HIERARCHY 

Let p k denote the result of the k th order relaxation (6) 
of problem (1), and let y k be the corresponding minimizer 
2fc t?I -order moment vector. 

Now, let us assume that the Archimedean condition (2) 
is met, and call r the degree of the polynomial on the 
right-hand side of eq. (2). Expressing the polynomials 

as MX) = Y. v fi v ( x ), 9ij ( x ) = T, v 9ijV( x ), 
follows that 


Cy \ tu = EE(/”)7^ t„w+ 

1=1 i v,w 

+EE(^)*^ tu w, (9) 

i,j v,w,s 

for |it| < k — [|"|. Due to positive semidefiniteness 
of the moment and localizing matrices, the right hand 
side of the above equation is non-negative, implying that 
Cy k tu > 2/«tx :fu’ f° r a ii Xl ' induction, it follows that 

> y k utu (10) 

for all sequences |u| < k — [|]. Such moments corre¬ 
sponds to the diagonal entries of the moment matrix 


Mk(y k ). Since Mk(y k ) > 0, it follows that ly^il < C^ w ^ 2 
for all monomials \w\<2k — 2\^~\. 

Now, for each vector y k replace with zeros all entries 
y k , with |ui| > 2k — 2|~I[], and complete the resulting 
2fc-order moment vector to an oo-order moment vector 
by placing even more zeros. We arrive at an infinite se¬ 
quence y s ,y s+1 ,... of vectors, with \y^\ < for all 

k,w. By the Banach-Alaoglu theorem [22], this sequence 
has a converging subsequence [70], and we will call y the 
corresponding limit. 

y satisfies in = 1, and Mk{y), Mk(qiy) > 0 for all k,i. 
By successive Cholesky decompositions of Mk(y), for k = 
s, s + 1,..., we find a sequence of complex vectors (|n))u 
with the property y u t„ = (it|u), for all monomials u,v. 
Call H = span{|u) : it}. We define the action of the 
operator X t on this (non-orthogonal) basis by 

Xi\u) = \X iU ). ( 11 ) 

and extend its definition to span{|tt) : u} by linearity. To 
prove that this definition is consistent, we need to show 
that, ifEu c “|n) = Eu^l^) f° r two different linear com¬ 
binations (c„) u , (d u ) u , then E U C «I^ M ) = E« d «l x iu). 
Indeed, note that, for any vector | w), 

(w| y^ y c u \Xju) = ^2,c u {w\Xiu) = ^CuiXiwlu) = 

u u u 

= (Xiw\^2c u \u) = {X i w\'^2,d u \u) = (wl^dulXiu), 

u u u 

(12) 

where the 2 nd and 5 th equalities follow from (w\Xiu) = 
y w \x iU = V(Xiw )tu = (Xiw\u). This relation holds for 
arbitrary \w), so the vectors E u c “|^n), J2 u ^\ X i u ) 
must be identical. Similarly, it can be verified immedi¬ 
ately that X.j is a symmetric operator, since {u\Xi\v) = 
Vu tx iV =y* v i Xi u = (' u\Xi\v) 

From the positive semidefiniteness of the localizing ma¬ 
trices M(qiy), it can be shown that {(p\qi{X)\(p) > 0 for 
all | (p) 6 H and i = 1, The Archimedean condition 

implies, moreover, that {(p\Xf\<p) < C(<p\<p). From this 
observation it is trivial to extend the action of Xi to 'H, 
the closure of H , and hence we arrive at a Hilbert space 
H and a set of operators X\, ...,X n such that qt{X) > 0 
for i = l,...,m and ij u = {ip\u(X)\'ip) , for |^) = |I). 

Note as well that, by construction, these operators sat¬ 
isfy all MPIs for dimension D. Now, call A the von 
Neumann algebra generated by Xi, ..., X n . By von Neu¬ 
mann’s 1949 result [23], such an algebra must decompose 
as a direct integral of type I, II and III factors [24] . That 
is, 

A = j dfj 1 (y)A l y © J dfi u {y)A l y ® j dfj lu {y)A l y l . 

(13) 
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Type I factors are isomorphic to B(H ), for Hilbert 
spaces TL of finite or infinite dimensionality [24]. Since 
A must satisfy the MPIs for dimension D, that excludes 
Hilbert spaces of dimension d > D from the first term of 
the right-hand side of (13). Moreover, in the Appendix it 
is proven that Type II and III factors violate the standard 
identity (5) for all values of d. 

It follows that we can write our operators Xi,...,X n 
as 

Xi = y* dfi I (y)Xi >y , (14) 

where each Xi tV acts on a Hilbert space H y with 
dim('H y ) < D. From qi(X) > 0, it follows that qi{X y ) > 
0 for i = 1 Hence p is a convex combination of 

feasible values of {p{X)) and so p < p*. On the other 
hand, p k > p* for k > s. Thus p = lim fc _ ) , 00 p fc > p*, 
proving the convergence of the hierarchy. 

Remark 2. Note that we just invoked the Archimedean 
condition (2) to establish the existence of y and, later, 
the boundedness of the operators X \,..., X n . Both results 
also follow from the weaker D-dimensional Archimedean 
condition; 

* 3 

+ Y J 9ij{X)Ui{X)g ij {X) + h D (X), (15) 

i,3 

where Hd{X) is an MPI for dimension D. 

Remark 3. If we take D = 1, then the MPIs will force 
all operators X\,...,X n to commute with each other. In 
that case, the SDP hierarchy reduces to the Lasserre- 
Parrilo hierarchy for polynomial minimization [25, 26]. 

Remark 4. So far, we have been assuming that the vari¬ 
ables Xi ,..., X n are Hermitian. If a subset of them is not, 
one can still define a converging SDP hierarchy, by con¬ 
sidering all possible monomials of the variables X\,..., X n 
and their adjoints x], ..., in the definition of moment 
vectors and moment matrices. In that case, the left-hand 
side of the D-dimensional Archimedean condition (15) 
must be replaced by C — JA 9QXJ — XjXi. 

V. EXPLOITING POLYNOMIAL 
CONSTRAINTS 

It is a basic result in operator algebras that MPIs for 
D x D matrices must have degree at least 2D [19]. This 
implies that we would need to implement the D th re¬ 
laxation of (6) or (7) in order to obtain non-trivial 19- 
dimensional constraints. Even for problems involving a 
small number of noncommuting variables, this becomes 


impractical already for D = 5. Hence, if we wish to con¬ 
duct optimizations over matrices of dimensions greater 
than 2 or 3, we must rely on linear restrictions other 
than those derived from MPIs. 

Most NPO problems relevant in quantum information 
science involve polynomial identities rather than polyno¬ 
mial inequalities. That is, constraints of the form q(X) > 
0 are complemented with —q(X) > 0, and so q(X) = 0 
must hold for all representations of {qi(X) > 0}. The 
strategy we will follow to solve this kind of problems is 
to divide the representations of Xi,...,X n into different 
classes r, in such a way that any two representations 
belonging to the same class r can be connected by a con¬ 
tinuous trajectory of feasible representations in r. As 
we will see, each of these classes will satisfy non-trivial 
low degree polynomial identities, which we will translate 
into linear constraints at the level of moment matrices 
(vectors). By carrying out a relaxation of the form (6) 
for each possible class r and taking the greatest result 
we hence obtain an upper bound on the solution of the 
general problem (1). 

For instance, suppose that {q^X) > 0} contains rela¬ 
tions of the form Xf = 1, for i = 1,2,3. For D = 2, there 
are two possibilities: 

1. Xi = ±1 for some i £ {1,2,3}. 

2. Xi ^ ±1 for all i, in which case it can be shown 
that the operators satisfy the identities: 


[X 1 ,{X 2 ,X 3 }+] = [X 2 ,{X 1 ,X 3 }+] = 

= [X 3 ,{A 1 ,X 2 }+]=0. (16) 

In either case, the noncommuting variables satisfy non¬ 
trivial polynomial constraints of degree smaller than 4, 
the smallest possible degree of an MPI for D = 2. A way 
to attack this problem is therefore to define an SDP re¬ 
laxation of the form (7) for each case, enforcing the corre¬ 
sponding extra linear constraints on the moment matrix 
(or moment vector). 

Note also that, if we further assume that the matrices 
X\,...,Xi are real, then we can add constraints of the 
sort 

{X 1 ,[A 2 ,X 3 ]}+ = 0 (17) 

in the second case. This approach hence allows us (in 
principle) to distinguish between real and complex ma¬ 
trix algebras. 

The objective, again, is to identify all possible linear 
restrictions on Mk or y within a given class r. Fortu¬ 
nately, most NPO problems in quantum information sci¬ 
ence have the peculiarity that random representations of 
a given class can be generated efficiently. 

Continuing with the previous example, suppose that 
we wish to optimize over 6 dichotomic operators, i.e., 
the polynomial constraints are, precisely, Xf = I, for 
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i = 1,2,...,6. For simplicity, let us denote the first 
four operators as X 0 o, X 0 i, Xi 0 , X 1:L and the last two by 
Yi , Yi . We want to maximize the average value of the 
operator 

p( x )=J2 E (—i rxc^YjXc^. (is) 

.7=1,2 ci,C2=0,l 

Note that we can write each dichotomic operator as 
Xj = (—l) a ~ g ° where {E l a } are projection operators 
satisfying E l 0 + E\ = I. Substituting in (18), we have 
that the objective function (^|p(X)|^) is equal to 

4 Y, -P( s > c i\X c1 c 2 jYi) + P(s, C2\X CiC2 ,Y 2 ) — 16, 

Cl,C2,S = 0,l 

(19) 

where P(ai, a 2 |xi, X 2 ) = E^\ijj). This cor¬ 

responds to the temporal correlations scenario defined 
in [18], where sequential dichotomic projective measure¬ 
ments are conducted over a quantum system, and a 
record of the measurements aq,X 2 ,... implemented, as 
well as the measurement outcomes asi, 02 ,is kept. The 
goal is to limit the statistics P(ai, ...,a n \x\, ...,x n ) ob¬ 
tained after several repetitions of the experiment. 

As noted in [14, 18], when the dimensionality of the 
quantum system is unrestricted, the set of all feasible dis¬ 
tributions P(ai, ...,a n \xi, ...,x n ), and hence the optimal 
value of (18), can be characterized by a single SDP. Us¬ 
ing the SDP solver Mosek [27], we find that, for D = 00 , 
p* = 8 up to 7 decimal places. 

Suppose, however, that we have the promise that the 
system has dimension D = 2. The problem we want to 
solve is therefore 

p* = max (ib\p(X)\ib) 

n,x,tp 

s.t. dim (%) < 2,1 — X? = 0, for i = 1,..., 6. (20) 

We start by dividing the representations of 2- 
dimensional dichotomic operators into classes. For any 
dichotomic operator, the rank of the projector E = 
can be 0, 1 or 2. For r = 0, 2, the corresponding operator 
is X = —I or X = I, respectively. For r = 1, a random di¬ 
chotomic operator X can be generated as X = 2 — I, 

where v 6 C 2 is a random complex vector. Since we 
are dealing with six non-commuting variables, there are 
3 6 = 729 classes, labeled by the vector r € {0,1,2} 6 , 
with rank(W + I) = Ti. 

For a fixed value of r, we sequentially generate random 
6-tuples of dichotomic operators X- 7 = (X(, ..., Xg), with 
the required rank constraints, as well as a sequence of 
random normalized vectors |i/d) £ C 2 . As before, we use 
each pair (X J , \tf>^)) to generate a random feasible mo¬ 
ment matrix Ml. Note that, since the conditions X 2 = I 
are implicit in each moment matrix, it is not necessary to 
include localizing matrices in our description (they would 


amount to zero diagonal blocks in the extended moment 
matrix). Notice as well that, given a feasible moment ma¬ 
trix Mfc, its complex conjugate M£ is also feasible. Since 
p{X) is a real linear combination of Hermitian monomi¬ 
als, the objective function will have the same value for 
both Mk and (and thus for the real feasible moment 
matrix Re(Mfc) = ^ This implies that, in order 

to define an SDP relaxation for (20), it suffices to con¬ 
sider the sequence of real matrices Re(M^), Re(M|),.... 

Applying the modified Gram-Schmidt method to that 
sequence until we find linear dependence, we obtain an 
orthonormal basis for A4 -, the space of all real feasi¬ 
ble moment matrices for representations of the class r. 
This time, the fact that this randomizing method works 
with probability one is a consequence that the projection 
of the randomly generated moment matrix Re(M^ +1 ) 
onto the orthogonal complement of the space spanned 
by Re(MjJ), ...,Re(M^) is a matrix whose entries are ra¬ 
tional functions of the randomly generated vectors used 
to build Xi +1 and \i/P +1 ). If Re(M^),..., Re(M^) do not 
span M k D p , then there exists a choice for those vectors 
such that the projected matrix is non-zero, i.e., at least 
one of such rational functions is non-zero. It is a well- 
known fact that the probability that a randomly evalu¬ 
ated non-zero rational function vanishes is zero. 

Alternatively, we can identify <Sp -, the space of fea¬ 
sible 2fc-order moment vectors for representations in the 
class r by parametrizing each normalized vector needed 
to build X or ip by two angles (f>, <p, and constructing the 
corresponding moment vector y(cj),0). Then the entries 
of the matrix 


S= j d<fdipRe(y($,ip))Re(y{<f,0)^) (21) 

can be computed analytically. Its support will coincide 
with Spf. 

One way or the other, we must solve the program 


p k =ma xY p w M w j, 

s.t. (22) 

for all possible classes r. For k = 2, again using 
Mosek [27], we find p 2 « 5.656854, definitely smaller than 
the free limit. 


VI. SIMILARLY INSPIRED SDP HIERARCHIES 

In the following we will introduce two problems in 
quantum information science which, while not exactly 
fitting in the class of problems (1), can be similarly re¬ 
duced to SDP hierarchies. 
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A. Quantum nonlocality under dimension 
constraints 

The scenario is as follows: two distant parties, call 
them Alice and Bob, conduct measurements on a bipar¬ 
tite quantum system. We denote Alice’s (Bob’s) mea¬ 
surement setting by x ( y ) and her (his) measurement 
outcome by a (6). We wish to bound a linear functional 
of the statistics P(a,b\x,y) they will observe, under the 
assumption that Alice’s and Bob’s spaces are, at most, 
.D-dimensional. If Alice and Bob’s outcomes are binary, 
i.e., a,b G {0,1} the problem can be shown equivalent to 

max B a,b P ( a ’ b \ x ’V)’ 

x,y,a,b 

s.t. P(a,b\x,y) = (ip\El <8> Fg\ip), (23) 

where {El,F^} are projection operators acting on C D , 
with £a El = F% = I D , and |^> G C D ® C D . 

Following the last section, we divide the representa¬ 
tions of the operators El , Fp into different classes labeled 
by the vectors r,t, with rank(D“) = r“, rank(D^) = t 
For each representation class r, t, we try to character¬ 
ize the span M . k ’ of feasible k th -order moment matri- 

D,rt 

ces. To do so, we sequentially generate random normal¬ 
ized states \ipi) G C D (g> C D and projectors G 

B(C d ), satisfying the rank conditions rank(D“’ ; ') = r“, 
rank(F b y ’ J ) = t%. Given El’i ,F^’\ we define the projec¬ 
tors = El'i^lo and F= I d®F^’\ which we will 
use to generate a feasible fc -order moment matrix M k . 
By subjecting the resulting sequence of moment matri¬ 
ces to the modified Gram-Schmidt orthogonalization, we 
obtain a basis for Ai k The SDP to solve is hence 

D,rt 

B k = max Y, B a,b( M k)E*,FZ 

x,y,a,b 

s.t. (Mfc)I,I = 1, M k > 0, M k G T k Pt (24) 

We again advise the reader to check that T = 
jr M J k is positive definite: otherwise, a projection 

of M k onto the support of T is necessary to guarantee the 
strict feasibility of the associated SDP, see Remark 1. 

The completeness of the above SDP hierarchy can 
be established easily: following the same lines as in 
Section IV, we prove the existence of a (in general, 
infinite-dimensional) representation El,F£ C B(P), 
with [E'l, Fff] = 0 for all x,y,a,b, and a state ip, such 
that J2x,v,a,b B a’,b (i’l^a^'bli’) coincides with the asymp¬ 
totic limit B = limfc_ s . 00 B k . The center of the algebra 
A generated by {El : x,a} decomposes W into a di¬ 
rect integral of sectors % z . By construction, in each sec¬ 
tor z, A boils down to a type I factor A z of dimension 
smaller than or equal to D. Being a type-I factor, we 
can write 1~L Z = Tif <g) Pf; then A z ~ B(Pp{) ® I, and 


A' z ~ I g) B(Pf), where A' z denotes the commutant of 
A z . It follows that 


K = j*drtz)EZ z ® Is,*, 


Fl = J® #)Ia, 2 ®^, 

(25) 


where dim (Pa,y) < D. Likewise, it can also be shown 
that the algebra generated by F^ z decomposes as a di¬ 
rect integral of finite dimensional algebras with dimen¬ 
sion smaller than or equal to D. B is thus a convex 
combination of feasible points and as such it represents 
a lower bound for the original problem (23). 

Examples 

Here we give examples of maximizing the violation of 
bipartite Bell inequalities with binary outcomes for differ¬ 
ent dimensionality of the component spaces. Some of the 
examples have already appeared in Ref. [15]. First we dis¬ 
cuss the 73322 inequality, the only tight three-setting two- 
outcome Bell inequality and its modified version. Then 
we move on to one more settings per party. For all the 
subsequent computations we used the solvers Mosek [27] 
and SeDuMi [28] through the interface YALMIP [29], 
which we ran on a memory enhanced desktop PC (with 
128 GB RAM). 

13322 — First we considered the 73322 inequality [30], 
which is the member of the Inn 22 family N > 2. Re¬ 
cently, it has been proven that qubit systems are not 
enough to attain the overall quantum maximum 0.2509. 
Rather, the best value in C 2 x C 2 systems is 0.25 [31, 32]. 
Using SDP, we reproduced the maximum value of 0.25 in 
dimensions C 3 x C 3 as well up to 8 significant digits [15]. 
The size of the moment matrix was 76 involving 1240 
linear constraints. The computations took about 5 min¬ 
utes for a fixed rank combination of measurements. Note 
that the hierarchy of Moroder et al. [32] , by limiting the 
negativity [33] of the bipartite quantum state, also gives 
a (not necessarily tight) upper bound on the Bell viola¬ 
tion for a fixed dimension of the quantum state. Indeed, 
this method works for C 2 x C 2 systems by returning a 
violation of 0.25 of the 73322 inequality [32]. However, for 
C 3 x C 3 systems it does not seem to converge (see Fig. 1 
of [32]). 

The SDP method also allows the user to upperbound 
the maximum quantum violation of a Bell inequality 
using a fixed two-qudit state. Let us choose the 3- 
dimensional maximally entangled state, \ip) = (100) + 

111) + |22))/x/3- We find the value of 0.229771, which is 
saturated by see-saw computation, hence the presented 
upper bound is tight. The computation involved a 116 
dimensional moment matrix (on a partial 4-level relax¬ 
ation) along with 1060 linear constraints. The program 
took 2 minutes to complete for a fixed rank combina¬ 
tion of projective measurements. We mention a related 
problem, where the maximal violation of 73322 has been 


computed for the maximally entangled state (of unre¬ 
stricted dimensionality). This has been solved both an¬ 
alytically [34] and using a relaxation method [35] by re¬ 
turning the value of 0.25. 

Modified 13322 — Though, I3322 inequality likely re¬ 
quires infinite dimensions to achieve the maximal vio¬ 
lation 0.2509 [13], D = 12 seems to be the smallest local 
dimension surpassing the qubit bound 0.25 [36]. It is 
an open question whether there exists a Bell inequality 
for which the maximum quantum violation in a given 
dimension is a strictly monotonic function of the dimen¬ 
sion. Below we give such a candidate. To this end, we 
modify the I3322 inequality by introducing a parameter 
c > 1: 

/3322(c) =Ef + E* + £f + £ 2 S 

— (£1, 1 + E\ 2 + £2,1 + £2,2) 

+ c(£i , 3 + £ 3,1 ~ -® 2 >3 — £ 3 , 2 ) < 4c, (26) 

where the correlator E XtV between measurement x by Al¬ 
ice and measurement y by Bob is defined as E x<y = P(a = 
b\x,y) — P(a ^ b\x,y), a,b £ {0,1} and Ef} denotes the 
marginal of Alice’s measurement setting x (and E y is 
similarly defined for Bob). This inequality is symmetric 
for exchange of Alice and Bob, and returns the original 
I3322 inequality (written in terms of correlators) for pa¬ 
rameter c = 1. 

Setting c = 2, we used the see-saw variational tech¬ 
nique [37, 38] to find a lower bound on the maximal vi¬ 
olation for any dimension 2 < d < 15, which we observe 
to be gradually increasing with dimension. We conjec¬ 
ture that the bounds are tight. Table I shows results up 
to D = 6 concerning both the lower (see-saw) and up¬ 
per bounds (SDP). Accordingly, the bounds for D = 2,3 
are indeed tight. Computationally the most challenging 
case was obtaining the upper bound in D = 4. It in¬ 
volved 3514 constraints, the dimension of the moment 
matrix is 184 and took roughly 40 minutes for Mosek to 
complete the task for a given rank combination of the 
measurements. The quantum maximum in dimensions 5 
and 6 coincide with the NPA bound on level 3 up to the 
shown digits. We pose it as a challenge to prove tight¬ 
ness of the see-saw bound for D = 4 (or possibly higher 
dimensions) by exploiting the symmetric structure of the 
inequality (26) using techniques such as in Refs. [32, 39]. 
14422 family - A one-parameter family of four-setting 
inequalities is given in Ref. [31]. These inequalities are 
not tight but they have a quite simple structure. They 
look as follows for c > 0: 

/4422(c) =c£j 4 + (/?i,i + E\2 + E 2 ,i — £2,2) 

+ (£ 43,3 + £ 3,4 + £ 4,3 — £ 4 , 4 ) <4 + c, (27) 

When c = 0, it is a direct sum of two CHSH inequalities, 
hence maximum violation is attained with qubit systems. 
However, by setting c > 0, it may serve as a dimen¬ 
sion witness. In particular, for c = 1 (the value used in 
Eq. (19) of [31]), its maximal violation in C 2 x C 2 systems 


D Lower bound Upper bound 


2 

8.013177 

8.013177 

3 

8.024050 

8.024050 

4 

8.032766 

8.071722 

5 

8.039579 

8.075937 

6 

8.056714 

8.075937 


TABLE I: Quantum bounds for different local dimensions on 
the violation of the 73322 ( 2 ) inequality computed using see¬ 
saw search/SDP computation. Bounds for D = 2,2, are tight, 
since the see-saw and SDP bounds match. As an overall upper 
bound, the NPA hierarchy on level 3 gives 8.075937. 


is upper bounded by the value of 5.8515 [31]. However, 
using our SDP tool, this upper bound turns out to be 
not tight: We certify a smaller value of 5.8310, which is 
matched by the see-saw method. Further, by raising the 
dimension to C 3 x C 3 , we get the same amount of viola¬ 
tion. The SDP computation returning 5.8310 in C 3 x C 3 
was quite demanding: it required a 130 dimensional mo¬ 
ment matrix and took about 2 hours of computational 
time. 5.8310 must be compared to the maximum value 
of 2v / 2 + v / 10 ~ 5.9907, achievable in C 4 x C 4 systems. In 
contrast to our certified value 5.8310, the corresponding 
C 3 x C 3 value arising from Moroder et al. hierarchy [32] 
(on their level 2) is a higher value of 5.9045. 

14122 inequality — It is also worth mentioning a situation 
(actually, this is the only case we are aware of) for which 
a previous SDP method introduced in Ref. [31] outper¬ 
forms our present SDP method. We tested the method 
in case of asymmetric Bell inequalities, that is, when the 
number of settings on the two sides are not the same. For 
the sake of comparison, we have chosen a correlation-type 
Bell inequality from [40], already analyzed in [31]: 

£4722 = £11 + £21 + £31 + £41 

+ (£12 - £22) + (£ 31 - £33) + (£41 - £44) 

+ (£ 25 - £35) + (£ 26 - £46) + (£37 - £47) < 8, 

(28) 

which consists of 4 and 7 binary-outcome settings on Al¬ 
ice and Bob’s respective sides. In Ref. [31] a method is 
presented in Sec. IIIB, which is particularly suited to 
asymmetric Bell setups. This way, the best upper bound 
obtained for qubit systems is 10.5102, whereas the best 
lower bound value of 10.4995 is due to see-saw search. 
The quantum maximum, attainable with two ququarts, 
is 10.5830. Using our present SDP technique and a desk¬ 
top PC, unfortunately, we did not manage to go below 
the global maximum 10.5830. 

Suppose now that Alice and Bob are conducting non¬ 
binary measurements, that is, measurements with more 
than just two outcomes. Then, in order to consider the 
most general measurements they could perform, we must 
model their measurement devices via Positive-Operator 
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Valued Measures (POVMs), rather than projective mea¬ 
surements. There are two ways to accomplish this: 

1. We can replace constraints of the sort ( E%) 2 = 
E% in (23) by the positive semidefinite constraints 
E% > 0, at the cost of having to add the corre¬ 
sponding localizing matrices to (24). Although con¬ 
verging, this method does not seem to behave well 
in our numerical experiments. 

2. Alternatively, we can exploit the fact that any 
d-outcome POVM {E a > 0} C B(C D ) can be 
realized in an extended Hilbert space C“ ® 
via a projective measurement of the form M a = 
C/(|a)(a| ® I d)U\ where U £ B(C d ® C D ) is a uni¬ 
tary matrix [41], Indeed, taking the state to be 
p = 10X01 ® IV'XV’I and choosing U appropriately, it 
can be verified that 


tr (pM a ) = tr(K|V , XV’l), (29) 

for a = 0,..., d — 1 and all states 

In the hierarchy to implement, random states of the 
form |0)(0|a' ® \i/j)(i/j\ab ® |0)(0|s' are generated. 
For each random state, we construct a moment 
matrix containing the operators E% = U x (\a)(a\ ® 
I D )(t/ x ) t ® I D ® Id, Fg = I d ® Id <8 V v (I D <g> 
\b)(b\)(V v y and the projectors Pa = |0)(0| ® I® 2 ® 
Id, Pb =Id®l| 2 ® |0)(0|. 

The convergence of this hierarchy follows from the 
fact that the algebras generated by {PaE^Pa}, 
{PbF^Pb} cannot violate D-dimensional MPIs. 


conjecture that the qubit bound 0.2071 using projective 
measurements is optimal (i.e. general POVM measure¬ 
ments do not help to improve the bound). However, in 
that case, we were unable to go below the known overall 
quantum maximum given in Refs. [13, 43]. 

The previous approach can be easily extended to char¬ 
acterize the statistics of multipartite scenarios where the 
local dimensionality of all parties is bounded from above. 
More interestingly, it can also be adapted to deal with 
multipartite Bell scenarios where only a subset of the 
parties has limited dimensionality. 

Consider, for instance, a tripartite scenario where Alice 
and Bob’s measurement devices are unconstrained, but 
the dimensionality of the third system (say, Charlie’s) 
is bounded by D. We want to generate a basis for the 
corresponding space of truncated moment matrices, with 
rows and columns labeled by strings of operators of the 
form u(AB)v(C), where u(AB) (v(C)) denotes a string 
of Alice and Bob’s (Charlie’s) operators of length at most 
kAB ( kc )• 

The key is to realize that, in a multipartite (complex) 
Hilbert space, the space of feasible moment matrices is 
spanned by moment matrices corresponding to separable 
states. Hence, in order to attack this problem, we start 
by generating a sequence of complex D-dimensional mo¬ 
ment matrices for Charlie’s system alone. After applying 
Gram-Schmidt to these complex matrices, we obtain the 
basis of Hermitian matrices {Mj}P =1 . Next, we generate 
a basis for Alice and Bob’s moment matrices. Since their 
dimension is unconstrained, such matrices are expressed 
as 

Fk AB = c uN u + CuN u t, (30) 

\u\<2k AB 


Examples 

We now apply our method to place nontrivial upper 
bounds on the quantum violation of Bell inequalities us¬ 
ing genuine POVM measurements for some of the set¬ 
tings. Note for binary-outcome settings general POVM 
measurements are not relevant, hence we have to con¬ 
sider Bell inequalities with at least one non-binary set¬ 
ting. To this end, we consider the simplest tight Bell 
inequality due to Pironio beyond genuine two-outcome 
inequalities [4, 42]. In this inequality, Alice has three 
binary-outcome measurements, and Bob has two settings: 
the first one has binary outcomes and the second one has 
ternary outcomes. If we allow Bob to use general POVM 
measurements on his second setting, the two-qubit quan¬ 
tum maximum (y/2 — l)/2 ss 0.2071 is recovered up 
to computer precision on level 3 of the SDP hierarchy. 
Hence, in this particular Bell inequality the use of gen¬ 
eral measurements do not provide any advantage over 
projective ones. Let us note that the quantum maximum 
without dimension constraints is a larger value 0.2532 
which can be obtained using a two-qutrit system and 
projective measurements [4]. We also applied the above 
method to the CGLMP inequality in order to prove the 


where N u is a matrix defined by 


(N u ) ViW = 1, if V s w = u ; 

0, otherwise. (31) 

The overall moment matrix for the whole system can then 
be expressed as M = ® ( c u,jN u + c* u j N ut ). 

Since we are just interested in optimizing a real linear 
combination of real entries of M —corresponding to the 
measured probabilities P(a,b, c\x,y, z) —, we can take 
the real part of the above matrix, and so we end up with 
the relaxation: 


max 

x,y,z,a,b,c 


B x f z M E x 

a,0,C a 




s.t. Mi,i = 1,M > 0, 

M = E ® ( N « + N u*)- 

U,j 

-c^.Im (Mj) ® (N u -N u t), 


(32) 
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where • (c^ •) denotes the real (imaginary) part of c u j. 
This is an SDP with real variables. 

Examples 

We now show applications of the above SDP method 
tailored to multipartite systems. As a first example, a 
three-party system is considered for which Alice possesses 
a qubit and the other two parties (Bob and Charlie) 
have no restriction on the dimensionality of the Hilbert 
spaces. We are able to fully reproduce the bounds ob¬ 
tained in Ref. [31]. In the next example, we extend Al¬ 
ice’s Hilbert space to a qutrit, thereby certifying genuine 
four-dimensional entanglement. Then we move to a four- 
party (translationally invariant) Bell scenario, and certify 
that a Bell value above a certain threshold cannot be ob¬ 
tained with symmetric measurements (that is, when each 
four parties measure the same observables in the first and 
second respective settings). 

1333 inequality — We consider the following three- 
party three-setting permutationally invariant Bell in¬ 
equality [31] 

/ 333 =sym{-P(A 1 ) - 2P(A 3 ) + P(A 1 ,B 1 ) 

- P(A 1 ,B 2 ) + P(A u B 3 ) - 2P(A 2 , B 2 ) 

+ 2P(A 2 ,B 3 )-2P(A 3i B 3 )} <0. (33) 

Here we used the short-hand notation P(A x ,B y ) = 
p(0,0\x,y), P(A X ) = p(0\x), and similarly for the other 
parties. Notice that the Bell expression above consists 
of only two-body correlators and single-party marginal 
terms, which usually provide an advantage in experi¬ 
ments. Such Bell inequalities have been proposed in 
Ref. [44] to detect nonlocality in multipartite quantum 
systems for any number of parties (however, those in¬ 
equalities involve only two settings per party, hence they 
can be maximally violated with qubit-systems, unlike the 
present example). In eq. (33), sym{A} means that ev¬ 
ery term occurring in X should be symmetrized with 
respect to all possible permutations of the parties, e.g. 
sym {P(Ai, Bi)} = P(A 1 ,B 1 ) + P(B U C{) + P(A U C{). 

We next compute upper bounds on the quantum vio¬ 
lations assuming different dimensionality of the Hilbert 
spaces. Lower bound values, on the other hand, are ob¬ 
tained from see-saw iteration in a prior work [31]. Ta¬ 
ble II summarizes the results. Values with an asterisk 
(*) have been established in the present work. Notation 
(DiD 2 D 3 ) refers to the dimensionalities of Alice, Bob 
and Charlie’s Hilbert spaces, respectively. Notice that 
due to symmetry of the Bell inequality (33), the same 
bounds apply to any permutations of (DiD 2 D 3 ). Estab¬ 
lishing upper bound on case (222) with non-degenerate 
measurements was the most time consuming task, the 
corresponding SDP problem involved 4894 constraints, 
and took 3 hours to be solved using Mosek—still the 
lower bound value has not been saturated. Comput¬ 
ing the upper bound for the case (2oooo) required to 
run the hybrid method (Alice was given level 2 of the 


LB UB 

LB 

UB 

LB 

(222) (222) 

(2oooo) 

(2oooo) 

(333) 

No-deg 0.0443484 0.0541362* 

0.1783946 

0.1783946* 

0.1962852 

Deg 0.1783946 0.1783946* 

0.1786897 

0.1786897 



TABLE II: Lower bounds (LB) and upper bounds (UB) on 
the violation of the I 333 inequality in various local dimen¬ 
sions. Values with an asterisk (*) have been established in 
the present work. The notation (D 1 D 2 D 3 ) refers to the di¬ 
mensionalities of Alice, Bob, and Charlie’s Hilbert spaces, 
respectively. The sign 00 denotes no restriction on dimen¬ 
sion of the respective party. Abbrevation Deg/No-deg refers 
to the situation when Alice has at least one degenerate mea¬ 
surement/all measurements are non-degenerate (i.e., rank-1 
projectors). The qutrit value (333) is the overall quantum 
maximum certified by the NPA hierarchy [12]. The upper 
bound value for (2oooo) in the degenerate case was obtained 
using the NPA hierarchy as well. 


qubit hierarchy, whereas Bob and Charlie’s system was 
computed on NPA level 1+AB). In that case, we man¬ 
aged to close the gap between the lower and upper bound 
values, thereby reproducing the result of Ref. [31]. We 
also computed (3oooo) upper bound and recovered the 
global maximum of 0.1962852 certified by the NPA hi¬ 
erarchy (Alice was given the level 3 of the qutrit hier¬ 
archy and took 24 hours for MOSEK to solve the re¬ 
sulting SDP). Accordingly, any Bell violation of I 333 big¬ 
ger than 0.1786897 cannot be attained with dimensional¬ 
ities ( 2 oooo) (plus the two other permutations) implying 
that the underlying three-party state pabc has at least 
Schmidt number vector (3,3,3) (see e.g. Refs. [45, 46]). 
Moreover any pure state decomposition of pabc contains 
at least one state uabc = such that the rank of 

each single party marginal < 7 ^ 4 , as, and ac is greater than 
2. In short, a Bell violation of I 333 bigger than 0.1786897 
detects in a device-independent way that the three-party 
state is genuinely three-dimensional entangled. 

1444 inequality — We construct a 3-party Bell inequality 
which cannot be violated maximally in state spaces C 3 x 
C D x C D (and arbitrary permutations thereof) for any 
dimension D. This extends the previous example to the 
case when Alice’s state space is restricted to a qutrit 
(instead of a qubit). In particular, the maximal violation 
is attained in C 4 x C 4 x C 4 . Hence, this certifies that 
the underlying three-party quantum state is genuinely 
four-dimensional entangled. 

Let us consider the following 3-party 4-setting Bell in¬ 
equality [47] 

/ 44 4 = CHSHab + CHSH^'C + CHSH B >c> < 6 (34) 

where A and A! denote different sets of measurements 
for party A, and we use similar notation for parties B 
and C. In Ref. [47] it has been proved that the maximum 
quantum violation attainable with biseparable states is 
S = 4 + 2\[2 ss 6.8284. Hence, if the above bound is 
exceeded in a Bell experiment, we can conclude that the 
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Bisep (2oooo) (3oooo) (444) 

6.828427 7.656854 7.971284 8.485281 

TABLE III: Maximum quantum bounds on Bell inequality 
7444 in (34) for different local dimensions of Alice. The first 
column (labeled by Bisep) stands for the case when Alice has 
a classical system and the other two parties have unrestricted 
dimensionalities. All bounds are tight as they are matched 
with lower bounds arising from see-saw iteration. 


state is genuinely tripartite entangled [48, 49]. The same 
bound can be derived by using the SDP techniques of 
Ref. [50] based on the NPA hierarchy. Below we extend 
this result to the realm of genuine higher-dimensional 
entanglement. 

To this end, we replace CHSH B <c/ with the Tsirelson 
bound 2v^2 [51]. This places an upper bound on 
I 444 in (34). Therefore, we are left with optimizing 
CHSH ab + CHSEU'c + 2^2 for C 3 x C D x C D sys¬ 
tems, where D denotes arbitrary dimension. To do so, we 
classify Alice’s four observables according to their traces 
(±1, ±3) and in each case we can solve the problem with 
SDP for the hybrid multipartite case. Notice, however, 
that Bob and Charlie have only two binary-outcome mea¬ 
surements, hence Jordan’s Lemma applies and we can as¬ 
sume that Bob and Charlie have traceless qubit observ¬ 
ables [52] . Then the problem goes back to upper bound¬ 
ing CHSH y 4B+CHSH J 4/c+2v / 2 in C 3 xC 2 xC 2 , which can 
be straightforwardly done using our SDP tools. By run¬ 
ning the SDP, the maximum turns out to be 36/7 + 2\/2 
up to the numerical precision of the solver Mosek. We 
also solved the problem assuming that Alice has a qubit 
yielding the upper bound 2 + 4#2 up to computer pre¬ 
cision. Results are summarized in table III. All bounds 
are tight as they are saturated using see-saw search. 

Consider now the so-called fully connected Bell state, 
that is, a 3-party state for which any two parties share a 
2-qubit Bell pair, 

1 ^ 444 } = \f + )ab ® \v + )a'c ® W + )b'C' ■ (35) 

With this particular C 4 xC 4 x C 4 state and measurement 
settings optimal for CHSH violation, we get the overall 
quantum maximum of 6>/2 « 8.485281 for the Bell in¬ 
equality (34). By adding a certain amount of white noise 
to the state (35): 

Pnoisy = P |^ 444 ) ('04441 H ^4x4x4? (36) 

we get the critical visibility p C rit = (36/7+2-\/2)/(6-\/2) ~ 
0.939425, above which we can detect the state (36) to be 
genuinely four-dimensional entangled. We believe this 
threshold is low enough to be interesting from an exper¬ 
imental point of view as well. 

12222 with symmetric measurements — In Ref. [53], a 
search has been conducted for all three- and four-partite 


binary-outcome Bell inequalities involving two-body cor¬ 
relators that obey translationally symmetry. Any trans- 
lationally invariant Bell inequality is provably maximally 
violated by a translationally invariant state when all par¬ 
ties measure the same set of observables (of unlimited di¬ 
mensionality). Numerical investigations in Ref. [53] sug¬ 
gest that it is not true anymore if we restrict the local 
Hilbert space dimension of the parties. Let us pick #64 
inequality from table 2 in Ref. [53] . Due to the fact that 
these Bell inequalities involve two dichotomic measure¬ 
ments per site, Jordan’s lemma applies and the maximum 
violation is given by (3q = 6 + 2v / 2 in qubit systems. Due 
to numerics, this value is achieved with different pairs of 
qubit observables. Indeed, running our SDP program by 
building up the bases from random symmetric measure¬ 
ments, we certify that applying the same settings at all 
sites do not allow us to violate the Bell inequality #64 
in table 2 of Ref. [53] (i.e., (3 q 1 = /3 C in the notation of 
the corresponding reference). 


B. 1-way quantum communication complexity 

Consider the following communication scenario: Alice 
and Bob are given inputs x, y with probability p(x, y ) and 
have the task to compute the Boolean function f(x,y). 
To do so, we allow Alice to transmit a D-dimensional 
quantum system to Bob, who, upon receiving it, must 
make a guess b on f(x,y). We wish to find the strategy 
which will allow Alice and Bob to maximize the probabil¬ 
ity that Bob’s guess is correct, i.e., b = f(x, y). For exam¬ 
ple, in a Quantum Random Access Code (QRAC) [54], 
the inputs x,y can take values in {0, l} fe and {1, 
respectively, and the function to compute is f(x, y) = x y . 

This scenario can be modeled by assuming that Alice 
prepares a pure quantum state p x = |# r XV’x| G B(C d ) 
depending on her input x. Bob will conduct a two- 
outcome projective measurement labeled by y, and de¬ 
fined by the projection operators {F]f : b = 0,1}, whose 
outcome will be Bob’s guess. In sum, we need to solve 
the problem 

max E p(x,y)tr(p x F^ x ^), 

s.t. tr(p x ) = 1 ,pl = Px , (if) 2 = if, 

p x ,Fy&B(C D ). (37) 

This problem can be reformulated by assuming that 
the initial state of Alice’s system corresponds to p x -o 
and, for any other input x, she sends the state V x poVJ. 
where V x is a unitary operator that can be chosen self- 
adjoint, i.e., V x = H. The resulting problem belongs to 
the class (1), and hence there is a converging SDP hier¬ 
archy to attack it. We observed that, in practice, such 
an SDP hierarchy gave good predictions for 77 = 2 at 
k = 2. For D = 3, a third-order relaxation did not suf- 






12 


fice to reach the optimal probability of success in 2 —> 1 
Quantum Random Access Codes (QRAC) [54]. 

We believe that the main reason for such a slow conver¬ 
gence rate is that the above proposal relies solely on MPIs 
to enhance dimension constraints. In order to devise 
a practical SDP hierarchy for problem (37), one needs 
to find a reformulation of problem (37) where the space 
Ai k D - is dimension-dependent even for low values of k. 
One such reformulation is immediate: regard p x as rank- 
1 projectors and assume that the state of the system is 
the (not normalized) tracial state I o- 

The resulting hierarchy of SDPs should be easy to 
guess: first, we divide the representations of problem (37) 
into different classes r depending on the rank of the pro¬ 
jectors {Fg}. Second, we generate random states p x and 
projectors {F//} within the class r, and, taking the state 
of the system to be Id, we use them to build random 
feasible moment matrices. Those allow us to character¬ 
ize the space A4/j -. Note that dimension constraints on 
A4 k D f are present for all D even for k = 1. E.g.: for any 
feasible first-order moment matrix M, Mjj = D x M\ tPx . 

The above SDP hierarchy gives good results in prac¬ 
tice, but we were not able to prove its convergence. Fol¬ 
lowing Section IV, it can be shown that, for any class r, 
one can define a representation F/f, p x and a tracial state 
which recover the limiting value of the hierar¬ 
chy of SDPs. Furthermore, the operator algebra decom¬ 
poses into a direct integral of representations z with the 
property that rank(.F, y ) < r]j, rank(I — p x ,z) < D — 1, 
rank(p XiZ ) < 1. If the dimensionality of H z is D , that 
defines a feasible point of problem (37). However, if the 
dimensionality of T~L Z is strictly smaller, we run into trou¬ 
ble: in such representations, p XjZ can vanish for some val¬ 
ues of x. Constraints such as tr (p z ) = 1 can be accounted 
for by other representations t where p X} t is a rank-1 pro¬ 
jector, since D{//) t \p Xit \%j t ) = > 1. 

One possibility to suppress the effects of lower finite¬ 
dimensional representations is to add ‘noncommuting 
constants’, i.e., certain extra operators whose operator 
relations cannot be realized in dimensions lower than D. 
For instance, in order to guarantee that all representa¬ 
tions Y have dimension D = 2, we could include the 
Pauli matrices a z , a x as operators in the moment matrix 
Mfc. With these extra variables, proving convergence can 
be done by appealing to the convergence of the Lasserre- 
Parrilo hierarchy [25, 26]. However, we did not find a sin¬ 
gle situation in our numerical experiments where adding 
noncommuting constants to fix the dimension was of any 
advantage. 

Examples 

We explore how the relaxation of the above communi¬ 
cation problem performs in practice. To do so, we estab¬ 
lish (usually tight) upper bounds in QRAC for various 
values of k and dimension D. We also recompute quan¬ 
tum bounds for the witnesses In of Gallego et al. [55]. 
We further distinguish between real and complex Hilbert 


spaces and detect general POVM measurements assum¬ 
ing that Alice communicates Bob a quantum system of 
fixed dimension D = 2. Note that Ref. [56] investigates a 
generalized QRAC problem where Alice’s inputs x take 
values from a string of dits (instead of bit-strings). In 
that case, our SDP method also showed good perfor¬ 
mance [56]. 

QRAC — We suppose the QRAC has independently 
and uniformly distributed inputs and Alice is allowed 
to transmit Bob a D-level quantum system. We use 
the notation of Ref. [57] and we denote the average suc¬ 
cess probability of the optimal k log 2 (D) QRAC by 

f’max(fc -» log 2 (£>)). 

It was known previously from Ref. [54] that P ma x( 2 —> 
1) = 1/2 + x/2/4. This is actually the value given by our 
SDP code at order 2, up to numerical precision. Likewise, 
when Alice is allowed to transmit a qutrit (case D = 3), 
our relaxation based on tracial states at the same order 
2 gives Pmax < 0.90450850, which matches with high nu¬ 
merical precision the lower bound value obtained via see¬ 
saw technique. Another method from the literature to at¬ 
tack this problem is the Mironowicz-Li-Pawlowski (MLP) 
SDP hierarchy [58], whose second-order relaxation gives 
us the (non-tight) upper bound of 0.9268355. 

One can prove that the MLP hierarchy does not con¬ 
verge in general. To do that, first notice that any QRAC 
can be rewritten as a full-correlation Bell inequality, by 
defining Alice’s observables as A x = 2 p x — I. Then, the 
only constraint that the MLP hierarchy adds to the NPA 
hierarchy is that (A x ) =2 —D. Taking D = 2, we see that 
the problem of calculating P max (fc —► 1) reduces to max¬ 
imizing the violation of a full-correlation Bell inequality 
constraining (some of) its marginals to be uniform. By 
Tsirelson’s theorem, the maximum is anyway attained 
when the marginals are uniform, so the constraint is au¬ 
tomatically satisfied [51, 59]. This means we can simply 
solve Tsirelson’s SDP to find out this maximum and the 
minimal dimension necessary to attain it [59]. For k = 4, 
we see that D — 4 is necessary to reach the maximum, 
so the hierarchy did not converge to the maximum for 
D = 2. 

By increasing the dimension D and the parameter k, 
the 2nd-order relaxation of the hierarchy based on the 
tracial states also performs well. The entries in the first 
three rows of table IV are from Ref. [15], which shows 
lower and upper bounds on the average success probabil¬ 
ity for QRAC k —> log 2 (D) for k = 3 and for different val¬ 
ues of D. The upper bounds (UB) are computed via our 
SDP in a normal desktop, and took less than 1 hour for 
any of the D values (assuming a given rank-combination 
of measurements) using the solver SeDuMi [28] . The up¬ 
per bounds (UB’) are resulting from the second-order 
relaxation of the MLP method [58]. We also show up¬ 
per bounds (UB”) derived from the Moroder et al. [32] 
hierarchy by fixing negativity (D — l)/2 and adding the 
constraint P(a\x) = 1/D on Alice’s marginal distribu¬ 
tions. As the table IV shows, except for D = 2,4, where 
the outputs of all methods coincide, the new tool gives 
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D 2_3_4_5_6_7 

LB 0.788675 0.832273 0.908248 0.924431 0.951184 0.969841 

UB 0.788675 0.832273 0.908248 0.924445 0.954123 0.969841 
UB’ 0.788675 0.853553 0.908248 0.934264 0.957785 0.979567 
UB” 0.788675 0.852156 0.908248 0.931201 0.954140 0.977072 

TABLE IV: Lower (LB) and various upper bounds (UB, UB’, 
UB”) on P max (3 —± log 2 (P)) detailed in the text. 


predictions ~ 10 -2 more accurate than the MLP method 
and the method based on Moroder et al. hierarchy. 

Let us pick P max (3 —> 1) = 0.788675 from table IV. 
This is precisely the value given by the construction of 
Ike Chuang [54] proving optimality of the complex qubit 
value. However, we can apply in this case the same ideas 
to characterize the properties of real qubit systems as 
well. By generating the basis from randomly chosen real¬ 
valued qubit states p x and projectors {Fq }, we get the 
(tight) upper bound 0.7696723. Hence, this simple ex¬ 
ample allows us to distinguish between real and complex 
two-level systems. 

In family — As another example, we used the SDP pro¬ 
gram based on tracial states to re-conrpute the maximal 
quantum value of the prepare-and-measure dimension 
witnesses In defined in Ref. [55], table I. The second re¬ 
laxation of the SDP hierarchy turns out to produce upper 
bounds for cases N = 3,4 and D = 2,3 which match the 
lower bounds obtained with see-saw method. Let us note 
that the conclusions of the experimental paper [60] relied 
on the conjecture that the inequality I 4 < 7.9689 cannot 
be violated by quantum systems of dimension D = 3. 

In a recent experimental paper [61], the In dimen¬ 
sion witness has been investigated for N = 7. Using 
a heuristic search, lower bounds are provided for dimen¬ 
sions D = 2,..., 6 , which were conjectured to be optimal. 
First row (LB) in our table V shows lower bound results 
of Ref. [61], which match our lower bounds (LB’ in sec¬ 
ond row) except the case D = 3, where we got a slightly 
higher lower bound value. This also justifies the need 
to relaxation methods providing certified upper bound 
values. Remarkably, our SDP upper bound values (UB 
in third row) are close to the LB’ values differing in the 
worst case D = 3 in the 2nd digit. For instance, Ref. [61] 
provides the experimental value of I 7 = 25.44 ± 0.02 for 
the preparation of 6 dimensional states. Therefore our 
UB value of 24.8991 for D = 5 certifies the generation of 
at least 6 dimensional quantum states. 

POVM witnesses — N. Gisin [62] asked if there exists a 
Bell inequality which requires POVMs for optimal viola¬ 
tion on some quantum state. This question has been an¬ 
swered affirmatively in case of two-qubit states (see e.g., 
Refs. [16, 17, 63, 64]). However, the question is still open 
for Bell inequalities defining facet of the local polytope 
(though, numerical study suggests the existence of such 
cases for three parties and high dimensional states [65]). 
We pose a similar question in the prepare-and-measure 


D 2_3_4_5_6 

LB 17.3976 20.7143 23.2167 24.8978 26.1017 
LB’ 17.3976 20.7085 23.2167 24.8987 26.1017 
UB 17.3976 20.7718 23.2180 24.8991 26.1019 

TABLE V: Lower and upper bounds on I 7 witness for dimen¬ 
sions D = 2,... , 6 . LB stands for the data in Ref. [61], table 
I, whereas LB’ is due to our see-saw technique and UB is 
resulting from our SDP method. 


communication scenario: By fixing dimension (say, Al¬ 
ice is allowed to send Bob a two-level system), are there 
witnesses which allow higher violations when Bob per¬ 
forms general POVM measurements instead of standard 
projective measurements? Our SDP tools allow one to 
certify the existence of such POVM witnesses. 

To this end, we pick the 14 witness of Ref. [66] and 
consider dimension D = 2. This witness consists of four 
preparations (x = 1,..., 4) on Alice’s side and six mea¬ 
surements (y = 1,..., 6) on Bob’s side. For D = 2, the 
maximum value of 14 equals 2-\/6, which can be attained 
if Alice prepares four states pointing toward the vertices 
of the regular tetrahedron (corresponding to SIC POVM 
elements). We add a four-outcome measurement to Bob 
(y = 7) with four outcomes b = 1,2,3,4 to the original 
14 witness and define the modified 14 witness as follows 

4 

Vf = 14 -J2 P ( b = A x = i,y = 7)< ~ 4.8990. 

i =1 

. (38) 

We remark that a similar modification was used in the 
context of Bell nonlocality in Ref. [63] . As the last term 
in the inequality cannot be positive, the qubit bound 
2^/6 ~ 4.8990 using POVMs follows from the bound on 
14. Indeed, by using the known optimal qubit settings 
for 14, the only way of getting the maximal violation of 
2^/6 is when Bob’s POVM elements in setting y = 7 are 
anti-aligned with the four tetrahedron states prepared by 
Alice, so that all probabilities P{b = i\x = i,y = 7) be¬ 
comes zero. By assuming projective qubit measurements 
for Bob and running SDP in case of D = 2, we obtain 
2(y/2 + 1) ~ 4.8284 up to numerical precision on the wit¬ 
ness Vf in Eq. (38). Hence, any value bigger than 4.8284 
for Vf certifies in a semi-device-independent way that 
Bob’s measurement y — 7 was in fact a general POVM 
measurement. 


VII. SOME TIPS ON IMPLEMENTATION 

In this section, we offer some tips to implement the 
programs defined above. As we will see, those tricks will 
make the hierarchy for NPO under dimension constraints 
much easier to code and modify than its dimension-free 
counterpart [14]. For simplicity, we will assume that our 
program does not involve localizing matrices, i.e., that 
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all polynomial restrictions appear as identities. We will 
also presume that operator representations can be di¬ 
vided into classes r, and that generating a random in¬ 
stance of each class can be done efficiently. 

Firstly, we need a subroutine 

[X,rho]=genSamp(D,n,r) that generates a cell-array 
X of random operators Xo, X\ : ..., X n , with Xq = Id 
and Xi, X n £ B(C D ) satisfying the appropriate class 
constraints, determined by the vector r. genSamp must 
also return a quantum state p £ B(C D ) (in our exam¬ 
ples, either a pure random state \ip)(ip\, the maximally 
entangled state or the unnormalized maximally mixed 
state). 

Secondly, we need a subroutine G = buildG(X,rho,k) 
that, given the cell-array X and an index fc, generates a 
k th -order moment matrix of the form: 

G u = tr(pXj k ...X}X jl ...X jk ), (39) 

where i,j £ {0, ...,n} fe , if the variables Xi,...,X n are 
Hermitian, or i,j £ {0,..., 2 n} k if they are not. In ei¬ 
ther case, monomials of G can be accessed via tr(G|i)(j|), 
choosing i,j appropriately. Note that, in this represen¬ 
tation, different columns of G correspond to the same 
operator, e.g.: 01 and 10. That does not matter, be¬ 
cause, at the end of the day, redundant columns will be 
suppressed by the matrix V mapping the space where G 
is defined to the support of span{G}, see Remark 1. 

If the reader is an Octave or MATLAB user, the fol¬ 
lowing considerations will lead to very fast code. 

Any D x D matrix A = i Ay I*) O'I can be repre¬ 
sented in vector form as |A) = J2i,j Aj^ji)01 : to go from 
A to |A), one can invoke the in-built function reshape. 
It can be verified that (Id ® {ip + \ ® Id)|A)|D) = |AD), 
where | i/j + ) is the non-normalized maximally entangled 
state | ij} + ) = J^iLi I M)- Similarly, (Id®( 0 |)|A) = 
and, hence, (A|(Id ® P*)\B) = tr (A^Bp), for any Hermi- 
tian matrix p £ B(C D ). 

Now, for random operators X \,..., X rl in the class r, 
define A = |Xj)(f|. From the above, it follows that 


(Id ® (ip + \ ®Id) I A® E \ x h- x n)(h,-ji\\ 

\ 32 ) 

= E \ X n- X h)(ju-M- (40) 

iu—Ji 

By tensoring A sequentially and projecting on the max¬ 
imally entangled state, we thus (quickly) obtain the op¬ 
erator 

G= E \ X n-X jk )(h,...i k \. (41) 

3l,---Jk 

Then it can be verified that 


G = G 1 '(Id ® P*)C. (42) 

Finally, we must code a third subroutine 
[basisMat,V]=buildBasis(D,n,r,k) that, by re¬ 
peatedly calling genSamp and buildG, derives an 
orthonormal matrix basis for span{G} and the isom¬ 
etry V described in Remark 1. From there, it is 
straightforward to implement program (7). 

Notice that all operator relations are determined by 
genSamp only. This means, for example, that if we 
wish to optimize over projection (unitary operators) all 
we need to do is program genSamp to generate ran¬ 
dom tuples of projection (unitary) operators P\,...,P n 
(Ui ,..., U n , U\ ,..., t/4). Switching from one type of poly¬ 
nomial constraints to another can thus be done straight¬ 
forwardly with the above implementation. 


VIII. CONCLUSION 

In this paper we have extended the notion of noncom- 
mutative polynomial optimization (NPO) to scenarios 
where the dimensionality of the spaces where the non¬ 
commuting variables act is bounded from above. We 
have presented a complete hierarchy of SDP relaxations 
to solve such problems and we have explored its perfor¬ 
mance by applying it to solve a number of open problems 
in quantum information theory. 

Our research raises several questions which deserve fur¬ 
ther study. The first one is whether the hierarchy of 
relaxations proposed to study 1-way quantum complex¬ 
ity is complete or not. As we showed, the main obsta¬ 
cle to prove convergence lies in interference from low¬ 
dimensional operator representations with ‘dark states’. 
For D = 2, the extreme points of such additional one¬ 
dimensional representations are finite, and thus the prob¬ 
lem of determining whether the hierarchy converges in 
a given prepare-and-measure scenario amounts to prov¬ 
ing that all such points can be reproduced by two-level 
quantum systems. Similarly, the convergence of the se¬ 
quence of relaxations to the maximum value of a specific 
functional for arbitrary D can be verified by establishing 
an upper bound for the dark-state value smaller than or 
equal to the value of a concrete D-dimensional realiza¬ 
tion. It would be more satisfactory, though, to have a 
general convergence result. 

Another open problem is whether the relaxations pro¬ 
posed for optimizations over finite-dimensional real oper¬ 
ator algebras actually converge. Here we are faced with 
the problem that certain operator representations, irre¬ 
ducible in the real space, can be expressed as a direct sum 
of non-trivial representations if we allow complex unitary 
transformations. Hence, if we wished to follow the proof 
for complex algebras, we would encounter problems at 
the step of applying von Neumann’s direct-integral the¬ 
orem [23]. 
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Finally, note that in this work we have not exploited 
the symmetry of the functionals to optimize. If our aim 
is to bound quantum nonlocality under dimension con¬ 
straints, that leaves us with a method that, in the best 
scenario, would allow us to conduct optimizations for 
D = 2,3 with a normal computer. Which dimension¬ 
alities could become accessible if we chose to play with 
symmetries is an intriguing question. 
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Appendix A: Type II and III factors violate the 
standard MPI 

The purpose of this Appendix is to prove its title. To 
do so, we will require the following lemma: 

Lemma 1. Let A be a type II or III factor. Then, for 
all n, there exists non-trivial projectors {Pi}" =1 C A such 
that 

1 • E"=l P i = L 

2. PiAP 2 A...P n A ± 0. 

Proof. Suppose that the statement is not true for general 
n and let TV be the greatest number such that it holds. 


Note that TV > 1. Indeed, if P\AP 2 = 0, for non-trivial 
Pi , P 2 then, for all x £ A, 


[x, Pi] = (Pi + P 2 )[x, Pi] (Pi + P 2 ) = 

= PixPi + P 2 xPi — PixPi — PixP 2 = 0. (Al) 

That is impossible, because factors, by definition are cen¬ 
tral, i.e., the only elements of A commuting will A are 
multiples of the identity. 

Now, suppose that {Pi }£Li satisfy the conditions of 
the lemma. Then we can always write Py = P' N + 
Pjv+i, where P' N , P ^ +1 are non-zero projectors such that 
P^AP^A.-.P^A ^ 0. Call B the algebra generated by 
AP^AP^A.-.P^A, and let n denote its identity, i.e., a 
projector n £ B such that Ay = y for all y £ B. 
Note that, due to the definition of B, AB = BA = B, 
and hence xll, Ax £ B for all x £ A. It follows that 
Ax = n^n = *n, that is [x, n] = 0. 

On the other hand, TV is the greatest number 
such that the conditions of the lemma hold, and so 
AP[...AP' n AP' n+1 = P' N+[ AP[...AP' N A = 0. In other 
words: BP' N+1 = 0, and, consequently, nP^ +1 = 0. We 
conclude that n, which is neither 0 nor the identity, must 
commute with all x £ A, contradicting the centrality of 

A. 

□ 


Now, given n and a type II or III factor A, choose 
orthogonal projectors {Pi}" =1 C A and operators 
xi,...,x n £ A such that PiXiP 2 x 2 ...P n ^ 0. Then de¬ 
fine the operators Ei = P^P^i, for i = 1,..., n — 1 and 
compute the fundamental polynomial I n -i(E). Due to 
the orthogonality of {Pi}, the only non-vanishing product 
is EiE 2 ..,E n —i = PiXiP 2 x 2 ...P n ^ 0. A hence violates 
In —1 = 0. 
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